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Abstract 

We solve the gl(l|2) generalized model by means of the algebraic 
Bethe ansatz. The resulting eigenvalue of the transfer matrix and the 
Bethe ansatz equations depend on three complex functions, called the 
parameters of the generalized model. Specifying the parameters ap- 
propriately, we obtain the Bethe ansatz equations of the supersym- 
metric t-J model, the Hubbard model, or of Yang's model of electrons 
with delta interaction. This means that the Bethe ansatz equations of 
these (and many other) models can be obtained from a common alge- 
braic source, namely from the Yang-Baxter algebra generated by the 
gl(l 1 2) invariant 7?-matrix. 
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1 Introduction 

An important class of exactly solvable 1+1 dimensional quantum systems is re- 
lated to the Yang-Baxter algebra. The theory of these systems is called the quan- 
tum inverse scattering method [0. The central object of the quantum inverse 
scattering method is the monodromy matrix, whose elements generate the Yang- 
Baxter algebra, and whose trace, the transfer matrix, determines a lattice statis- 
tical model, which, in many cases, is of interest in physics. The structure of the 
Yang-Baxter algebra is fixed by a set of numerical functions, the elements of the 
7?-matrix, which satisfies the famous Yang-Baxter equation [EH. 

Solving an exactly solvable model means at first instance to calculate the spec- 
trum and the eigenfunctions of the transfer matrix. In many cases of interest this 
task can be accomplished using the quadratic commutation relations between the 
elements of the monodromy matrix. The model is then said to be solved by alge- 
braic Bethe ansatz. For the algebraic Bethe ansatz to work it has to be possible 
to identify the elements of the monodromy matrix as 'particle' creation and an- 
nihilation operators. In particular, a pseudo vacuum state must exist, which is 
annihilated by all the annihilation operators. In all known cases, where an alge- 
braic Bethe ansatz was successful so far, the elements of the monodromy matrix 
could be arranged in such a way, that the monodromy matrix acts as an upper 
triangular matrix on the pseudo vacuum, and the pseudo vacuum is an eigenstate 
of its diagonal elements. Still, even if a pseudo vacuum exists, on which the mon- 
odromy matrix acts triangularly, no general recipe for an algebraic Bethe ansatz 
is known, and the remaining calculations may be rather involved (see [^] for a 
highly non-trivial example). 

The algebraic Bethe ansatz is best understood for models with the 7?-matrix of 
the spin-^ XXX and XXZ chains [[]]]. In this case there is only one monodromy 
matrix element below the diagonal, which must annihilate the pseudo vacuum. 
The particular form of the vacuum eigenvalues of the diagonal elements of the 
monodromy matrix is not used in the construction of the transfer matrix eigenvec- 
tors. The vacuum eigenvalues, usually denoted a(X) and d(X), enter the algebraic 
Bethe ansatz solution as functional parameters. Thus, one may think of a model 
defined by the functional parameters and the triangular action of the monodromy 
matrix on the pseudo vacuum. The question whether such kind of 'generalized 
model' can be realized for arbitrary parameters a(k) and d(k) was first addressed 
in [f|] and after some refinement was answered affirmatively in For the 

models with XXX and XXZ 7?-matrix the arbitrariness of the functional parame- 
ters was the key tool in Korepin's proof [|7p of the norm formula and in Slavnov's 
work on form factors. 

The simplest models which allow for a nested algebraic Bethe ansatz are the 
models with gl(n) invariant 7?-matrix [R IM. Considering the fundamental rep- 
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resentations of these models one observes that not only the monodromy matrix 
elements below the diagonal annihilate the pseudo vacuum, but additional zeros 
appear above the diagonal [JIT]]. This fact simplifies the algebraic Bethe ansatz 
for the fundamental representation as compared to the more general case, where 
the action of all the elements of the monodromy matrix above the diagonal is 
non-trivial. For the solution of this more general case a new concept, the vacuum 



subspace, was introduced by Kulish and Reshetikhin Q12Q. This new concept en- 
abled to perform the algebraic Bethe ansatz for the models with gl(n) invariant 
7?-matrix with the same generality as in the gl(2) case. The resulting eigenvalue of 
the transfer matrix and the Bethe ansatz equations depend on n functional parame- 
ters ai(X),...,a n (X), which, together with the triangular action of the monodromy 
matrix on the pseudo vacuum, define the so-called gl(n) generalized model [fT3|] . 
Considering the parameters a\ (X) , a2(A,) , a^(X) as free parameters Reshetikhin de- 



rived the norm formula for the gl(3) generalized model [|T3|]. 



In this work we construct the algebraic Bethe ansatz solution for the gl(l|2) 
generalized model. This work has three motivations: 

(i) To perform the algebraic Bethe ansatz for all models with gl(l |2) invariant 
graded 7?-matrix [§, |14|] . 



(ii) To pave the ground for the calculation of the norm for this class of models. 

(iii) To point out an interesting relation to the Hubbard model, namely, that the 
Lieb-Wu equations JT3p can be generated by the Yang-Baxter algebra with 
gl(l 1 2) invariant 7?-matrix. 

Let us comment here on point (i) above. Points (ii) and (iii) will be further dis- 
cussed in section 7. 

Numerous articles have appeared about models connected to the Yang-Baxter 
algebra generated by the gl(l|2) invariant 7?-matrix. In first place we have to 
mention the lattice gas model of Lai and Sutherland [[16|,^7|] and the super sym- 
metric t-J model [ITS]] . Their relation to the graded Yang-Baxter algebra was 



explored in [|19|, [20|] . Both models are connected to the fundamental represen- 
tation of gl(l|2). Their algebraic Bethe ansatz solution was obtained in JRj ] (see 
also [jl9p). More recently there was a wave of interest in models related to higher 



dimensional representations of gl(l |2) 021|,[22D. A fermionic model related to the 
four-dimensional typical representation was introduced in [|Z3"1-|Z5|] and was named 
supersymmetric U model. Algebraic Bethe ansatz solutions of this model were ob- 
tained in [^-^Bp (see also p9|]). An impurity t-J model, where the fundamental 



representation at one lattice site is replaced by the typical four-dimensional repre- 
sentation, was considered in [|3"U|, [31]] . Finally, the latest example in our list is the 
family of 'doped Heisenberg chains' [|32p related to the atypical representations of 
gl(l|2). 
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The algebraic Bethe ansatz solutions of all the above mentioned models can 
be obtained as special cases of the solution presented in this work. In this sense, 
we claim that we performed the algebraic Bethe ansatz for all models with gl(l |2) 
invariant i?-matrix. We would also like to point out that our result applies to the 
Bethe ansatz for the quantum transfer matrix of the supersymmetric t-J model 
and can be used to prove the conjectures, presented in [ [331 , [34] ], about the quantum 
transfer matrix eigenvalue and the corresponding Bethe ansatz equations. 

Let us note that one of the Bethe ansatz solutions in our above list is quite 
different from the others and is not covered by our approach. Ramos and Martins' 
solution P%| ] is based on the Yang-Baxter algebra generated by the intertwiner of 
two typical four-dimensional representations of gl(l |2) and not on the fundamen- 
tal 7?-matrix. They obtain the same result for the eigenvalue of the transfer matrix 
of the supersymmetric U model and the same Bethe ansatz equations as in []26|], 
but their construction of the eigenvectors is quite different. It will be interesting 
to find out the relation between the two seemingly different sets of eigenvectors. 

The plan of this work is the following: In section 2 we make more precise what 
we mean by the gl(l |2) generalized model. In preparation of the algebraic Bethe 
ansatz we work out the block structure of the Yang-Baxter algebra in section 3. 
Section 4 is the core of this work. It is devoted to the algebraic Bethe ansatz. 
In section 5 we consider the fundamental graded representation, related to the 
supersymmetric t-J model. This section has been included to give an example 
of how our slightly abstract formalism relates back to physics. Following the 
example the reader will be able to work out by himself the cases he is interested in. 
In section 6 we point out a relation of our general result to the Lieb-Wu equations. 
Section 7 is devoted to a discussion of possible applications. 

In order not to overload this article with formulae we shall restrict ourselves 

to a fixed choice (H ) of the grading. The other two possible cases (see e.g. 

[ |T4"[ |T9p) are left to the reader. They present no more difficulty than the case 
considered in the text. 



2 The gl(l 1 2) generalized model 

The starting point of our considerations is the gl(l |2) invariant rational i?-matrix 
[Fl|l4l] with matrix elements 

^)=a(X)(-l)^)8^ + H^)5^ , (1) 

a,(3,y,8 = 1,2,3. For simplicity we restrict ourselves to the grading p(\) = 0, 
p(2) = p(3) = 1. The matrix solves the Yang-Baxter equation and is in- 
variant with respect to the action of the fundamental representation of gl(l |2). It 
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satisfies the compatibility condition [Q] 

^^ ) = ( _ 1)P (a)+p(P)+p(Y)+ P (S)^Y (X) _ (2) 

The complex valued functions a (A) and b(k) of the spectral parameter A G C are 
defined as 

A + ic A + ic 

They satisfy the equation a (A.) + fr(A) = 1. The additional complex parameter c 
will be called the coupling constant. We shall further need the matrix/? (A) defined 
by 

^) =*]£M • (4) 

The graded Yang-Baxter algebra with i?-matrix R(X) is the graded, associative 
algebra (with unit element) generated by the elements 1jj a (A) , a, (3 = 1 , 2, 3, of the 
so-called monodromy matrix modulo the relations 

R(X-p)(T(X)® s T(^) = (T(p)® s T(X))R(X-p) . (5) 

We assume that the elements of the monodromy matrix are of definite parity, 
7t(Tp a (A) ) = p(ct) + p(P) . The symbol ® s denotes the super tensor product associ- 
ated with the grading p(l) = 0, p(2) = p(3) = 1. For a definition see appendix A. 

The gl(l|2) generalized model is the set of all (linear) representations of the 
graded Yang-Baxter algebra (Q) with highest vector Q, defined by the fact that 
T (A) acts triangularly on £1, 

T 1 1 (A)^ = ai(A)a , <r 2 2 (A)C2 = a 2 (A)a , <r 3 3 (A)Q = a 3 (A)C2 , (6) 
Tp a (A)^ = , for a > (3 . (7) 

The complex valued functions ay (A), j = 1,2,3, are called the parameters of the 
generalized model. These parameters characterize the representation in a similar 
manner as the highest weight in a highest weight representation of a Lie algebra. 

Let us denote the representation space of a given representation of the general- 
ized model by H . It is clear from the quadratic commutation relations contained 
in the graded Yang-Baxter algebra (|5|) and from (|(J), (0) that we may assume that 
H is spanned by all vectors of the form 

4»(A 1 ,...,A^) = T«>(Ai)...^(A w )^ , (8) 
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where < p\, k = 1, . . . ,N. This assumption is at least sensible for a finite 
dimensional representation space H . 
Let us define the transfer matrix 

t(X) = (-iyW<r ( ?(X) = str(<r(X)) . (9) 

Since R(X) is invertible for generic values of X E C, we conclude from (0) and @ 
that 

[t(X),t( M )} = (10) 

for all generic X,/j E C. Thus t(X) and have a common system of eigenfunc- 
tions. In other words, the eigenvectors of t(X) are independent of the spectral 
parameter X. The task of the algebraic Bethe ansatz for the generalized model is 
to diagonalize t(X), i.e., to solve the eigenvalue problem 

t{X)<f> = A{X)<f> . (11) 

It is a remarkable fact that this task can be accomplished by solely using the graded 
Yang-Baxter algebra (|J) and the properties (§) and (|7]) of the highest vector Q.. 
In particular, it is not necessary to require that T^(X)Q. = as in case of the 
fundamental graded representation, which corresponds to the supersymmetric t-J 
model. 



3 The graded Yang-Baxter algebra 

The structure of the graded Yang-Baxter algebra (|J) becomes much clearer after 
rewriting it in block form. We introduce the shorthand notations 



B(X) = {B l (X),B 2 (X)) , C(X)=(^ 2 ^j 

Then we can write the 3x3 monodromy matrix ? (X) as 

(A{X) Bi(X) B 2 (X)\ f 
<I{X)=[C\X) D\(X) D\(X))=(pV B{X) 

\C 2 (X) D\(X) D\{X) ^ C{A) D[k) 



(12) 



(13) 



The defining relations of the graded Yang-Baxter algebra (Q) can be thought of as 
a 9 x 9 matrix equation. Let us denote the n x n unit matrix by /„. A similarity 
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transformation with the matrix 



X 



(h 



\ 




\ 



h) 



(14) 



which cyclically permutes the 5th, 6th and 7th row, transforms this 9x9 equation 
into 



/l 



\ / A®A A®B B®A B®B \ 

bl 2 al 2 A®C A®D -B®C -B®D 

al 2 bl 2 C®A C®B D®A D®B 

f) \-C(g)C -C®D D(g)C D&D ) 



( A® A A®B B®A B®B \ 

A®C A®D -B®C -B®D 

C®A C®B D®A D®B 

\-C®C -C®D D®C D®D J \ 



(\ 



bh ah 
ah bh 



\ 
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• (15) 



For the formula to fit on the line we suppressed the arguments and adopted the fol- 
lowing convention: A,...,D depend on X, and a bar means here that X is replaced 
by /j. Furthermore, a = a(X — /u) and b — b(k — /j). The 4x4 matrix 



(b 



\ 



\ 



b —a 
—a b 



(16) 



b — aj 

is of 'six-vertex form' . It satisfies the Yang-Baxter equation and is unitary, 

f(X)f(-X)=I 4 . 



(17) 



We would like to remark that the defining relations of the graded Yang-Baxter 
algebra of the gl(l|2) model, when written in block form (JT3|), resemble the cor- 
responding relations for the gl(l 1 1) model. 

Out of the 16 relations contained in (jT5J) we shall need the following 4 for the 
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algebraic Bethe ansatz, 

B(X) ®B(ji) = (BQi) ®B(X))f(X-^) , 

w^n ~r s b(u-X) , v B(fj)®A(X) 

A(X)®B(/j) = - f- . j . B(X)®A(jj) + - v ' 



D(X)®B(jj) 



a (ji — X) 
b(k-n) 



a (ji — X) 



B(X)®D(fi)-(B(ji)®D(X)) 



r(X-») 



a(X — /u) a(X — n) 

r{X-n)(D{X)®D{n)) = (D{n)®D{X))f{X-n) . 



(18) 
(19) 

(20) 
(21) 



Note that, by (|T%[), B(X) constitutes a representation of the Zamolodchikov alge- 
bra, and, by (pit), D(X) is a representation of the Yang-Baxter algebra of the gl(2) 
model. 



4 The algebraic Bethe ansatz 

Our goal is to calculate the eigenvectors of the transfer matrix t(X) = A(X) — 
tr(D(X)). In analogy with the gl(l|l) case we shall first of all consider the com- 
mutation relations of a multiple tensor product B(X\) ® ■ ■ -®B(Xn) with A(X) and 
tr(D(X) ) . These commutation relations can be obtained be iterating equations ( fl^ ) 
and (|2"0"|). We first present the result of the iteration and explain our notation below, 



N 



j=1 
7=1 



N 



N 



1 



A(X) <g)B(Xj) = ®B(Xj)A(X)H 

L j=i J j=i a \ K i K ) 

N 



Hi 



D(X) 



N 



7=1 



B(Xj) =(-l) N \h 



N 



a(Xj-X) \ = \a{X k -Xj) 

Hi 

B(Xj)})f(X)fl- 



7=1 



f = \a(X-Xj) 



N 



N 



-(-l) N Y,\h®B(X)® <g)B(X k ) \PoiSf\ 

7=1 



Hi 



l 



a(X-Xj) f = \a(Xj-X k ) 

Hi 



(22) 



(23) 



Here the operators B(X) in the multiple tensor products are multiplied in ascending 
order. We make use of the usual conventions for the embedding of linear operators 
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into tensor product spaces and define 

r j - lJ (X)=lf j - 2) ®r(\)®lf N - j) (24) 

for j = 2, ... ,N as an operator in (End(C 2 )) 8>A '. The operators S 7 _i appearing on 
the right hand side of (|22|) are then defined as 

Sj-i = ri ) 2(A,i-A,j)r 2j 3(A,2-A, J ')...rj_ij(A,y_i-A,y) , (25) 

for j = 2, . .. ,N. We further define So = id. Then the notation used in ( ]2"2| ) is 
explained. 

It remains to explain the operators Pq\, 5y_j and r(X) appearing at the right 
hand side of equation (|2~3|). These operators act in one more auxiliary space 'zero'. 



sfl^h^Sj^ , P 0l =P®lf N - l) , (26) 

where P is the permutation matrix on C 2 g) C 2 with matrix elements P^g = 5g8p. 
In order to define f(k) we introduce the gl(2) i?-matrix 

r(X) = Pf(X) (27) 

and the monodromy matrix of the corresponding fundamental inhomogeneous 
model, 

T(°XX)=r^-^N)...r { °}(l-l l ) . (28) 

We shall interpret this auxiliary monodromy matrix as a 2 x 2 matrix in space zero 
with entries acting on spaces 1, . . . ,N. We further introduce 

D { ° ] (X)=D(X)®lf N . (29) 
Eventually, T(X) is defined as 

f(X) =D { ^(X)T^(X) . (30) 

This expression, too, can be understood as a 2 x 2 matrix, the trace of which 
appears under the sum on the right hand side of equation (C 



Equations and (g3j) were obtained by the usual symmetry arguments of 
the algebraic Bethe ansatz. Their derivation solely relies on (P"8])-(|ZU|). As an 
immediate consequence of (jT|) we obtain the equation 



N r N \ 

(g)B(k k ) = h(k j )® <g)B{X k ) \Sj-i , (31) 

t=1 ^ k=\ J 



k=\ K k=\ 
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which was repeatedly used in the derivation of (|2~2"|), (|2~3"|). A mathematically rigor- 
ous proof of (j22l), (|2~3l), which is not too hard to do, may be obtained by induction 
over N. Some useful identities needed in the proof of the less trivial case of ( |2"3| ) 
are provided in appendix B. 

A commutation relation, appropriate for our purposes, of the transfer matrix 
t(X) with a multiple tensor product B(Xi) <g> ■ ■ ■ ®B(Xn) now easily follows from 
( |2~2"1 ) and (p3|). We have to take the trace of equation ( |23| ) in space zero and have 
to subtract the result from equation (Q). Using the fact that b(X)/a(X) = ic/X is 
an odd function of X we obtain 



N 



t(X) ®B(Xj) 



N 



N 



7=1 



2V 



My 



1 



^aiX-X^UaiXj-Xk) 

k+i 

i=1 a{k k -kj) 

k^j 



(32) 



Before proceeding with the Bethe ansatz calculation we have to gain more 



insight into the mathematical structure of the various expressions in equation (|32|): 
The vector B(X\) (g> • • -®B(Xn) may be thought of as a 2^ component row vector 
with entries acting on the representation space H . The components of this vector 
are 



®B{Xj) 

.7=1 



11,...,/JV 



•JVJ 



(33) 



i n =1,2. Equivalently, ® • ■ -®B{X N ) is a vector in (C 2 )®^ ® End(#). 

Suppose F lu -'' N e ^ for all z'i, . . . , = 1,2. Then F can be defined as a column 
vector with 2^ components F ll '- lN that are vectors in or, equivalently, F e 
(C 2 )®^®^. It follows that 



r w i 

[(g)B(X J )\F=B h (X l )... B iN (X N )F'^ e X 

7=1 



(34) 



Following Kulish and Reshetikhin JT2] ] let us define the 'vacuum subspace' 9{q C 



11 



H by the conditions 

A{X)<f> = ai(X)<!> , (35) 

C(X)<f> = , (36) 

for all 4> G 9{q. Clearly, 9(q is a linear subspace of 9{ . The following lemma 



holds [13]. 



Lemma 1. H§ is invariant under the action ofD(X). 

Proof. From the Yang-Baxter algebra written in block form ( |T5| ) we deduce the 
equations 

-b(l-/Li)B(X)®C(n) +a(X-fi)D(X) ®A(ji) = 

a(X-/u)A(iu)®D(X)-b(X-iLi)B(iLi)®C(X) , (37) 

r(X-/u)(D(X)®C(v)) = 

-a(X-n)C(n)®D(X) + b(X-n)D(n)®C(X) . (38) 



Let 4> G 9{q. Acting with fl37j ) and ( |38D on $ we obtain 



(A( A /)(8)D(X))4) = fli( A /)D(X)4» , (39) 
(C{n)®D{X))<S> = Q , (40) 

and the lemma is proven. □ 

Corollary 1. The space spanned by all linear combinations of vectors of the form 
D\(/ui) . ..D\(hm) Q is a linear subspace of Mq. 

Thanks to our lemma we can now proceed with the so-called second level 
Bethe ansatz. We introduce the notation 

T {X) -[c(o\X) D«»(k)) ' T{X) -{c(X) D(X)J ■ (41) 

Then, by the definition @ of T(X), 

A(X) =a(°\X)d\(X) +C {0) (X)D\(X) , 

B(X)=B (0 \X)D\(X)+D^(X)D\(X) , 

C(X) =A {0 \X)D 2 l (X) +C {0 \X)Dl(X) , 

D(X) = B(°\X)D\(X) +D ( ~°\X)Dl{X) . 
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Here we used the fact that 

[D^(X),T^l(X)}=0 , fara,P,Y,8=l,2. (43) 

Since both, D(X) and (X), are representations of the Yang-Baxter algebra with 
7?-matrix r(X) (see (|2"T|)), equation ( |43| ) allows us to conclude that 

f(X-^)(f{X)(Bt{ M )) = (f{ M )(Bf{X))f{X-fi) . (44) 

It follows that tr(f (X)) can be diagonalized by the algebraic Bethe ansatz if a 
pseudo vacuum exists on which T(X) acts triangularly. 

The operators A(X),. . . ,D(X) act on (C 2 )®^ In this space let us define 
the vector 

6=(J)®%o . (45) 

We shall see that £i is an appropriate pseudo vacuum for T(X). Let us denote the 
elements of the gl(2) standard basis by e«. Taking into account that 

,(°)m_ (ip-a)ej\-aej\ be s \ \ 

r ^ k) ~\ bejj -aej\ + (b-a)ejl) ' (46) 

when written as a matrix in auxiliary space, and using the definition ( |28| ) of the 
auxiliary monodromy matrix, we obtain 

N 

A^(X)Q m = Yl(b(X-Xj)-a(X-Xj)).(y N , (47) 



N 



D m w© w =(-i) w n«(*-»*)-©"" • < 48 > 

d°)(X)(i)^ = o . (49) 



Hence, using @, (|7]) and d42|), we conclude that 

A(X)Cl = (-lfa 2 (X)f\ a -^^Cl 

j=1 a{kj-k) 



(50) 



D(X)a = (-lfa 3 (X)Y[a(X-Xj)a , (51) 
C(X)£l = . (52) 



This means that T(X) acts triangularly on Cl. 
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Let us briefly recall the algebraic Bethe ansatz solution of the gl(2) generalized 
model (see e.g. [[]]]). Out of the 16 relations contained in (Q) we pick out the 
following three, 



B(X)B( M )=B( M )B(!) , 



(53) 
(54) 

(55) 



Moving A (A) and D(A) through a product B(ji{) . . .B(/j m ) by means of (p3|)-(p5|) 
leads to 



M 



a (a) 



l k=l 



M M i 

M r M 



1 



¥* 



l 



. M M 
M . M 



1 



l^k l^k 

We add the latter equations and use the fact that b(X)/a(X) is odd. Then 



(56) 



(57) 



M 



(A(A)+D(A)) 



I 

yfc= 1 



/t=l 



M 



M 



k=l 

M 



i 



M 

[J a{X-Hk) 



M j 

a{n k -X) 



+ £ B(A)[]B(W) 



1 



(58) 



Now £2 is a joint eigenvector of A(A) and £>(A) (see (|5(j|), (pTj)). Thus, acting with 
both sides of equation (|5%|) on £2 we see that 5(^i) • ..B(hm) Cl is an eigenvector 
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oftr(f (A)) =A(X)+D(X), 

tr(T(X))B( Ml )...B( MM )& = A(X)B( M i)...B( MM )& , (59) 
with eigenvalue 

7=1 

if the Bethe ansatz equations 

«3 (A**) /J a(Xj - fit) i=i a (fn - jUfc) 

are satisfied for = 1 , . . . , M. 

Since the vacuum subspace H§ C H is invariant under D(A), it follows that 
fl(jui) . . .5(a/m) ^ e (C 2 )® w (8) #b- In other words, BQii) . . .5(/iM) ^ is a column 
vector with 2 W rows having vectors vo.Hq<ZH as entries. Hence, by the definition 
([35]), ( 06| ) of the vacuum subspace, B{ji\) . . .B(/jm) ^ is an eigenvector of A (A), 
viewed as an operator on (C 2 )®^ ® , 

A(A) J B( w )... J B( m )^ = a 1 (A) J B( A i 1 )...5( m )^ . (62) 

Recall that fi(Ai) <8> ■ ■ ■ ®B(Kn) is a row vector with 2 W columns. It follows from 
(02), © and © that 



(g)5(A y ) B{ni)...B{p M )& 

7=1 



= B/, (AO . . .B iN (X N ) [B(m) . . .B(p M )£l] iu ~' iN e X (63) 
is an eigenvector of t(X) with eigenvalue 



A(A)=ai(A)n 



n ! 



N y M l M 1 

M a(Ay - A) ^ a( A - n k ) ^ a - A) 
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if the Bethe ansatz equations 

are satisfied for j = 1, . . . , N. Thus we have obtained the eigenvectors (O) and 



eigenvalues ( |64| ) of the transfer matrix t(X) of the gl(l |2) generalized model. The 
eigenvalues and eigenvectors depend on two sets of Bethe ansatz roots {kj}j =1 , 
{rtj^Li determined by the Bethe ansatz equations ( |o"T| ) and (|53|). We would like 
to stress that we did not specify the functions aj(X), j = 1,2,3, yet, nor did we 
require that D^ityQ. = 0. Our calculation solely relied on the graded Yang-Baxter 
algebra (Q) and on equations @ and (|7]). 

Remark. Suppose all Xj and /j^ be mutually distinct, and the parameters aj(X) 
have no singularities at Xj, Then A(X), equation (|53|), has simple poles at 
X = Xj, j = 1, . . . ,N and X = /Jk, k = 1, . . . ,M. The condition for the poles at 
X = Xj to vanish is 

res{A(^)}=0 , j=l,...,N (66) 
and is equivalent to (j65l) . Similarly, the equations 

res{A( A /,)}=0 , fc=l,...,M (67) 

are equivalent to d6l)). 



5 The fundamental representation 



In order to illustrate our results with an example let us reconsider the algebraic 
Bethe ansatz of the supersymmetric t-J model [jT9p. We have to connect the model 
to a representation of the graded Yang-Baxter algebra (|]) and have to identify 
the parameters a\(X), ai{X) and a-s(X) of the representation. The Bethe ansatz 
equations are then given by (|6"T|), (|65|), the eigenvectors and the eigenvalues of the 
transfer matrix by ( |63| ) and (|64|). We shall basically follow the account of [ |35l , |36| ]. 

First of all, we introduce canonically anticommuting creation and annihilation 
operators cj a , Ck,b of electrons (a,b =f, j; j,k= 1, . . . ,L) and the Fock vacuum 
|0), annihilated by all annihilation operators, c^O) = 0. The elements (Xj)^, 
a, P = 1 , 2, 3, of the matrix 



X, 



(68) 



7i c ;t n i\i l - n n)J 
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form a complete set of projection operators on the space of states locally spanned 
by the basis vectors |0) , c1j|0) , cj| |0) . Double occupancy of lattice sites is forbid- 
den on this space. Let — (-Xy)g, oc, (3 = 1,2,3. The operator Xy™ = 1 — n^nj^ 
projects the local space of lattice electrons onto the space from which double oc- 
cupancy is excluded. The corresponding global projection operator is 



Po 



L 

rid 



n n n ji) 



(69) 



Owing to the fact that the operators are local projection operators, it fol- 
lows from general considerations [E51] that the 'L-matrix' 



fa 



Lj(X)=a(X)J 3 + b(X) 



v.l 
A J2 
V 2 v 2 
A Jl ~ A J2 



1 

2 
1 

Ji 



v.3 

A J2 



_v.2 
-^3/ 



(70) 



is a representation of the graded Yang-Baxter algebra @. This representation has 
been termed fundamental graded representation in [|35|]. The action of Lj(X) on 
the Fock vacuum obviously is 



Lj{X)\0) 



f\ b{X)Xj\ b(l)Xj\\ 
a(X) 
v y 



|0> 



(71) 



The matrix Lj(X) generates the supersymmetric t-J model at a single site. The 
corresponding monodromy matrix of the L-site model is 



<T(X)=L L (X)...L 1 (X) 

Its action on the Fock vacuum follows from ( [71] ) as 

(\ Bi(k) B 2 (X)\ 
a L (X) 



(72) 



<T(X)\0) 



|0) 



(73) 



\0 a L {X) J 



From the latter equation we can read off the parameters of the representation, 
ai(X) = 1, ci2(X) = a^iX) = a L (X). And we are done (compare JT^], equations 
(3.47), (3.48) and (3.50)). 



Remark. Note that the Hamiltonian of the supersymmetric t-J model is 

1 



H = -ic^lnj (str(T (0))) str(T (X)) } 



(74) 
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Because it acts on the restricted space of electronic states, where no lattice site is 



doubly occupied, we may replace it with (see Q35|, p6[|) 

L L 



HPo = Po{ - £ (c$ a c j+ha + cj +ha c ha ) + 2 £ (S?S* +l 



7=1 



4 



+ n 



Po 



(75) 



6 Connection to Lieb-Wu equations 

We demonstrate that for an appropriate choice of the parameters of the gl(l|2) 
generalized model the Bethe ansatz equations (|6T|), ( |65| ) turn into the Lieb-Wu 
equations [03J], which are the Bethe ansatz equations of the Hubbard model. 

Let us introduce two functions of the spectral parameter, k(k) and v(X), de- 
fined by 

s'mk = 'k , v = X — ic/2 . (76) 

Setting 

kj = k{Xj) , j=l,...,N , (77) 
v* = v(wt) , k=l,...,M (78) 

and using the explicit form (Q) of the function a(A,), we can rewrite the Bethe 
ansatz equations (^5J), ( |6T1 ) as 

ai(smkj) _ -it y fc - sin kj - ic/2 
a2(sinfcj) t^Vk — sinkj + ic/2 

a2(vk + ic/2) X Vfe-shifcj-ic/2 = v* - v/ - ic/2 
a 3 (vjt + ic/ 2) A 1 Vfe - sinkj + ic/ 2 }J v fc - v/ + ic/ 2 

for j = 1, . . . , /V; & = 1, . . . ,M. Obviously, the latter two equations turn into the 
Lieb-Wu equations through the choice 

«i(sinfc) jkL a 2 (v + ic/2) 
a2(sinfc) a3{v + 1C/2) 

of the parameters of the generalized model. 

Similarly, replacing sink with k in equations (f75|)-(|8"TD we obtain the Bethe 
ansatz equations [ f37| ] of Yang's model of electrons interacting via a delta function 
potential. 
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7 Conclusions 

We have obtained the algebraic Bethe ansatz solution of the gl(l|2) generalized 
model. The Bethe ansatz equations (|61|), (|65|) and the eigenvalue of the transfer 
matrix (|4|) depend on three functional parameters, which are determined by the 
respective representation of the graded Yang-Baxter algebra (|5|). Choosing the 
parameters appropriately we obtain the Bethe ansatz equations for the models dis- 
cussed in the introduction. In order to calculate the transfer matrix eigenvalues 
for the models based on higher representations of gl(l|2) we have to supply ad- 
ditional arguments based on the ideas of fusion and analyticity of the eigenvalue 
(see©). 

What are the further applications of our result? They relate to points (ii) and 
(iii) of the list in the introduction. We do not have much doubt that Reshetikhin's 
work QT3p on the norm of the gl(3) generalized model carries over almost liter- 
ally to the gl(l|2) case (see point (ii) in the introduction). It is, however, point 
(iii) which enthralls us most. We propose to systematically seek for a represen- 
tation of the graded Yang-Baxter algebra (|5|) with parameters a\(X), «2(^) an d 



03 (A,) satisfying (8jJ). Such a representation would certainly much improve our 



understanding of the algebraic structure of the Hubbard model. It would probably 



enable to prove the norm conjecture proposed in [ pBQ for the Bethe ansatz wave 
functions of the Hubbard model, and, combined with the results of [ffijj], might 
eventually lead to form factor formulae for local operators in the finite Hubbard 
chain. 

Unfortunately, for the case of gl(l|2) no theorem is known so far that would 
guarantee the existence of a representation of the graded Yang-Baxter algebra 

with arbitrary functional parameters. It would be highly desirable to find out 
whether or not a generalization of the work of Korepin and Tarasov to the 
gl(l|2) case is possible. As long as this question is not settled it remains unclear 
whether the appearance of the Lieb-Wu equations in the context of the gl(l|2) 
generalized model reflects a deep connection with the Hubbard model or is just a 
coincidence. 



Acknowledgments 



The author likes to thank V. E. Korepin and M. Shiroishi for helpful discussions. 
He thanks K. Sakai for encouraging him to publish his results. 



19 



Appendix A: Graded algebras 

In this appendix we shall recall the basic concepts of graded vector spaces and 
graded associative algebras. In the context of the quantum inverse scattering 
method these concepts were first utilized by Kulish and Sklyanin [|9|,[T4|]. 

Graded vector spaces are vector spaces equipped with a notion of odd and 
even, that allows us to treat fermions within the formalism of the quantum inverse 
scattering method (see [ f?5| , |3£i| ]). Let us consider a finite dimensional vector space 
V, which is the direct sum of two subspaces, V = Vq © V\ , dim Vq = m, dim V\=n. 
We shall call vo G Vo even an d v\ G V\ odd. The subspaces Vq and V\ are called 
the homogeneous components of V . The parity n is a function V{ — > Z2 defined on 
the homogeneous components of V, 

n( Vi ) = i , i = 0,l , Vi eV . (A.l) 

The vector space V endowed with this structure is called a graded vector space or 
super space. 

Let PL be an associative algebra (with unity), which is graded as a vector space. 
Suppose X,Y G PL are homogeneous. If the product XY is homogeneous with 
parity 

n(XY) =k(X)+k(Y) , (A.2) 

then PL is called a graded associative algebra [p. 

For any two homogeneous elements X, Y G PL let us define the super-bracket 

[X,Y}± =XY-(-l) n ^ n ^YX , (A.3) 

and let us extend this definition linearly in both of its arguments to all elements 
of PL. 

Let p : {1 , . . . ,n} — > Z2. The set of all nx n matrices A, 5, C, . . . with entries 
in PL, such that 71 (A?) = % (B?) = 7t(Cp ) = • • • = p(a) + /j((3) is an associative 

algebra, say Mat(j*,n), since 7t(Ap5y) = p(a) + p(y). For A, 5 G Mat(^,n) we 
define the super tensor product (or graded tensor product) 

(A® a fi)g=(-l)W°)+^W>WAj5j . (A.4) 

This definition has an interesting consequence. Let A,5,C,D G Mat(j?,«), such 
that [B« Cj]± = 0. Then 

(A®,5)(C®,D) =AC®,5D . (A.5) 
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Appendix B: Three identities 

In this appendix we provide three identities that are useful for the proof of equation 
Cp3b. The first one is 



N 



N >v 

(g) B(Xjj\ ® h \ {lf N ®D(X) } P N _ m . . . P l2 Poi 



■ Ax-x N )..4kx-x 1 ) 



N 



®B(Xjj\ \d^(X)T^(X) . 

7=1 



(B.l) 



The second one is 

. N 



Hi 



B(X)® <g)B(X k ) ®D(Xj)\r^l lN (X J -X N )...rf ) j+l (X J -X j+l ) 



*(0) 



7 2 ® [® | ^01 ■ • .P N -w{l? N ®D(Xj)} 

■ *n- l ,/v ( X i ~ • • ■ rfj + ! (fy - V i ) 



TV 



h®B{X)® (g)B(X k ) \P 0l ...P j - lj 

■Df\Xj)rfl(X j -X N )...rf) +1 (Xj-X j+1 ] 



= \l 2 ®B(X)® ®B(X k ) kfj^^-Vi)' 

Hi 

. ..rfliXj-XjD^iXjyfUXj-Xv) . ..rfl^Xj-Xj+r) 

( N } 

= ll 2 ®B(X)® [(g)B(A*)] ^i^l^^tro^^rC )^))} . (B.2) 



M7 
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The third identity is an identity between rational functions that can be easily 
proven by means of Liouville's theorem, 

b(k — X^+i) yj 1 
a(X-X N+ \) JLi a(k N+ i -X k ) 

_ b(X-X N+l ) » 1 N b(Xj-X N+ i)b(X-Xj) N 1 

a(X-X N+i ) M a(X-Xj) A a(Xj-X N+ i) a(X-Xj) l = \a(Xj-X k ) 

(B.3) 

The induction step in the proof of (23|) may be done as follows. First, insert ( |B.1| ) 
and QB.2P into (|23|). Second, take the tensor product of the resulting equation with 
B(Xn + i). Third, use ( |T8| ) and ( |20| ) to bring the operators into the required order. 



Finally, use (JO) to rearrange the resulting terms in the desired way. 
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